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Abstract. We study a matrix- valued Schrodinger operator with random point 
interactions. We prove the absence of absolutely continuous spectrum for this 
operator by proving that away from a discrete set its Lyapunov exponents do 
not vanish. For this we use a criterion by Gol'dsheid and Margulis and we 
prove the Zariski denseness, in the symplectic group, of the group generated 
by the transfer matrices. Then we prove estimates on the transfer matrices 
which lead to the Holder continuity of the Lyapunov exponents. After proving 
the existence of the integrated density of states of the operator, we also prove 
its Holder continuity by proving a Thouless formula which links the integrated 
density of states to the sum of the positive Lyapunov exponents. 



1. Introduction : model and results 

This paper discusses properties of the Lyapunov exponents and the of integrated 
density of states of the following formal Schrodinger operator with random point 
interactions, 



dx 2 



( ciu[ n) 5o(x-n) \ 

\ c N uj^S (x-n) J 



acting on L 2 (M.) ® C . Here, Vq is the maximal multiplication operator by the 
constant coefficient matrix with on its diagonal, 1 on the upper and lower di- 
agonals and everywhere else. Moreover, c±, . . . , c/v are real numbers, Sq is the 
Dirac distribution at the point and In is the identity matrix of order N, N > 1. 
Let (f2, A, P) be a complete probability space on which we define the sequence 

(u>^)nez — ((w| n ', ■ ■ ■ jk'jy ))nez of independent and identically distributed (i.i.d.) 
K^-valued random variables of common distribution v. We assume that the sup- 
port of ^, supp v C M. , is bounded and that {x — y | x, y G supp v\ spans R^. An 
example of such a distribution is obtained when the components , . . . , wjy ' are 
i.i.d. real- valued random variables such that the support of their common distribu- 
tion contains at least two different points (like in the case of Bernoulli variables). 
We also set (fi, A, P) = ®„ eZ (0, A, P). 

Following PQ, we define the self-adjoint realization of the formal operator Hp(u>), 
for every u> £ fl, by 

JV 

(2) HH=0fl w< +V o 

i=i 
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acting on L 2 (R) ® C . Each H Ui for i E {1, . . . , N} is the operator acting on 
L 2 (K) (g) C by HuJ = -f" whose domain is 

D(H^) = {/ g L 2 (R) ® C | /,/' are absolutely continuous on E \ Z, 

/" g L 2 (R) (g) C, / is continuous on E, 

f'(n + ) = f '(ri~) + CjCj| f(n) for every n g Z}, 

where existence of the left and right derivatives f'(n~) and f'(n + ) at all integers 
is assumed. The Wj's, Vo and the c%'a are as above. These operators are self-adjoint 
and bounded from below (see [I]) as the support of the distribution v is bounded. 
Thus, Vo being bounded and self-adjoint, it implies that H(u>) is self-adjoint for 
every uj g Q. Moreover, H(uj) is a Z-ergodic operator and therefore its almost-sure 
spectrum is included in E. 

The operator H(u>) is the Hamiltonian which describe the propagation of an electron 
in N coupled one-dimensional layers. The random variables are describing point 
sources that interact with the electron. These point sources are situated at each 
integer point of the system. 

In the case N — 2, a previous work from Giinter Stolz and the author (see [5]) 
stated that H(u>) has almost-surely no absolutely continuous spectrum. We now 
look at the case N = 3 which is our first result on model J5J). 

Theorem 1. We assume that N = 3 and that ci g E, c\ and C2 are non-zero. 
Then there exists a discrete set cScl such that for every E £ E \ S, the Lyapunov 
exponents associated to H(lo) at E are all distinct and positive, 

7 i(£) > 12 {E) > 73 (£) > 0. 

Thus, H(u>) has almost-surely no absolutely continuous spectrum. 

The approach to prove this theorem will be to study the Fiirstenberg group G(E) 
associated to every real number E to H(v). That is, the group generated by the 
transfer matrices associated to H(u>). We will recall precise definitions and proper- 
ties of the transfer matrices, the Fiirstenberg group and the Lyapunov exponents 
in Section [2~T1 Let Sp N (M) denote the group of 2N x 2N symplectic matrices. The 
Fiirstenberg group G(E) is a subgroup of Sp N (E) . For the definitions of L p -strong 
irreducibility and p-contractivity we refer to [2], definitions A/V.3.3 and A.IV.1.1, 
respectively. Then the proof of theorem [1] is based upon the following criterion due 
to Gol'dsheid and Margulis. 

Theorem 2 ( 10 ). If the group G(E) is Zariski dense in Sp N (E) then it is p- 
contracting and L p -strongly irreducible for every p g {1, ... , N}. Thus the Lya- 
punov exponents associated to H(tu) are distinct and positive. 

Once we have proved the Zariski denseness of the Fiirstenberg group, by Kotani's 
theory (see [14]) we will obtain the almost-surely absence of absolutely continuous 
spectrum of H{w). Actually [14] does not cover directly the case of point interac- 
tions models as ([2]) but we will see in Section 13.2.21 how to adapt proofs of [14] to 
our setting. 

Our second result will be on the regularity of the Lyapunov exponents viewed as 
functions of the real parameter E. We first prove a general result of Holder con- 
tinuity of the Lyapunov exponents associated to H(oj) under suitable assumptions 
on G(E). 
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Theorem 3. Let I be a compact interval in R. We assume that the group G(E) is 
p-contracting and L p -strongly irreducible for every p € {1, . . . , N} and all E 6 I . 
Then the Lyapunov exponents associated to H(uj) are Holder continuous on I , i.c, 
there exist two real numbers a > and < C < +oo such that 

Vpe{l,...AT}, VE,E'eI, \ 7p (E)- lp (E')\<C\E-E'\ a . 

The proof of theorem [3] is mostly based upon the existence of an integral rep- 
resentation of the Lyapunov exponents under the assumptions of p-contractivity 
and Lp-strong irreducibility of G(E). It will also require estimates on the transfer 
matrices that will be proved in Section 13.11 Combining theorem [T] and theorem [3] 
we also obtain, 

Theorem 4. We assume that N = 3 and that C2 € R, c\ and c-i are non-zero. Let 
S be the discrete subset of R obtained in theorem [7J Let I C R \ S be a compact 
interval. Then the Lyapunov exponents associated to H(uj) are Holder continuous 
on I . 

Finally in Section [321 we turn to the study of the integrated density of states. It 
is a function of one real variable that counts the mean number per unit volume 
of spectral values of H(uS) below a fixed real number E. For operators like H(ui), 
acting on an infinite dimensional Hubert space, this function could be equal to 
infinity for every real number E. To avoid this problem, we will first restrict the 
operator to intervals of R of finite length. Let L be a strictly positive integer and 
D = [-L, L] C R. We set H^ D \u) the restriction of H{w) acting on L 2 {D) <g> C N 
with Dirichlet boundary conditions on D. 

Definition 1. The integrated density of states, or IDS, associated to H(u>) is the 
function N : R — > R + defined for each £ £ I as the following thermodynamical 
limit, 

(3) N(E)= lim i#{A<£|A6 ff (ff< D )(w))}. 

In Section 13.2.11 we will start by proving the existence of this limit. For this, we 
will have to prove a matrix-valued Feynman-Kac formula for matrix-valued point 
interactions models. To prove the existence of the IDS for every E 6 R, we will 
adapt ideas of [7] for scalar-valued point interactions models combined with what 
was done in for matrix-valued Anderson models. Then in Section ^. 2. 2l we deduce 
a Thouless formula by adapting Kotani's theory to point interactions models. This 
formula states that the sum of positive Lyapunov exponents and the integrated 
density of states are harmonically conjugated. Combining this Thouless formula, 
properties of the Hilbert transform and theorem [3l we prove the following result. 

Theorem 5. Let I be a compact interval in R. We assume that the group G(E) is 
p-contracting and L p -strongly irreducible for every p £ {1, . . . , N} and all E G / . 
Then the integrated density of states associated to H(u>) is Holder continuous on I. 

Combining theorem [5] with theorem [1] we get this last theorem. 

Theorem 6. We assume that N = 3 and that C2 S R, c\ and C3 are non-zero. Let S 
be the discrete subset o/R obtained in theorem]^ Let I c R\<S be a compact interval. 
Then the integrated density of states associated to H(u>) is Holder continuous on I. 
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The fact that we can choose C2 = in theorem [T] and theorem [6j means that we 
are here in the presence of a phenomenon of propagation of randomness similar 
to the one found in the work of Glaffig (see [9]). Let us explain what we observe 
here. If C2 = or, as we will see in Section 12.21 if H^ 2 is deterministic, we still 
have the positivity of all the Lyapunov exponents and the regularity of the IDS of 
H(uS), But, due to Kotani's theory, the positivity of the Lyapunov exponents is 
directly related to the randomness in the model. Heuristically, if the second layer 
in our model (corresponding to H U2 ) is deterministic, we should not have that all 
Lyapunov exponents are positive. But the fact that we have randomness on the first 
and the third layer and that the three layers are coupled by the action of Vq seems 
to mean that the randomness is somehow transported to the second layer. To the 
best of our knowledge the only other example of such a phenomenon of propagation 
of randomness can be found in [9J where Glaffig proves a strong regularity result 
for the IDS of a discrete matrix-valued Schrodinger operator. In Glaffig's model, 
the randomness on the first layer is assumed to be very strong as it is given by i.i.d. 
random variables whose common law is very regular (absolutely continuous with 
respect to the Lebesgue measure and with a Radon-Nikodym derivative being in 
a Sobolev space). On the other layers he only assumed that the random variables 
follow a Bernoulli law. Despite that lack of regularity on the other layers, he 
still proves that the IDS is C°° which should have occured intuitively only if all 
the random variables were as regular as those on the first layer. But, as all the 
layers are coupled in his model, the strong regularity on the first layer is somehow 
propagating to the other layers. 

Even if we do not prove that the IDS in our model is C°° like in Glaffig's model, 
we prove the Holder continuity with even less randomness than in Glaffig's model. 
We only need Bernoulli variables on the first and the third layers and we can even 
assume that the second layer is deterministic. Still, the Holder continuity of the 
IDS is sufficient to hope to be able to prove in the future the Anderson localization 
for -ff(w), that is that the spectrum of H(ui) is pure point and the associated 
eigenfunctions are exponentially decaying to at infinity. Indeed, with the Holder 
continuity of the IDS we should be able to prove a Wegner estimate which is the 
key ingredient in order to apply a multi-scale analysis scheme to prove Anderson 
localization (see [T61I8]). 

We finish this introduction by giving the outline of the rest of the paper. Section^ 
will be devoted to the proof of theorem[TJ In Section [2~T1 we will recall the definitions 
of the Lyapunov exponents, the transfer matrices and the Furstenberg group. Then 
in Section \2. 21 we will present the computations which leads to the Zariski denseness 
of the Furstenberg group in the symplectic group. These computations will be 
very similar to those done in [TUj for the one-dimensional discrete matrix-valued 
Schrodinger operator. The reason is that in the point interactions model that we 
study here, the random parameters only appear at integer points. Thus, in the 
transfer matrices, the random parameters will appear in the same way as in the 
discrete model. The difference will be that we have here a continuous model and 
thus the energy parameter E will appear in a way different than in the discrete 
case. This is the reason of the existence of the set S of critical energies in theorem 
[I] These critical energies did not appear in the discrete model of [ID] . 
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In Section [3] we will focus on the regularity results for the Lyapunov exponents and 
the IDS. Section [3. II deals with the regularity of the Lyapunov exponents while in 
Section 13.21 we will prove the existence of the IDS and its Holder continuity. 

2. POSITIVITY OF THE LYAPUNOV EXPONENTS 

2.1. Lyapunov exponents and transfer matrices. We start this section with a 
review of some results about Lyapunov exponents. These results holds for general 
sequences of i.i.d. random symplectic matrices. Let N be a positive integer and 
Sp N (R) denote the group of 2N x 2N real symplectic matrices. It is the subgroup 
of GL2n(K) of matrices M satisfying t MJM = J, where J is the matrix of order 
-In ' 



2JV defined by J : 



In 



Definition 2. Let (A£) n ett be a sequence of i.i.d. random matrices in Sp N (R) 
with E(log + \\Aq II) < +oo. The Lyapunov exponents 71,.. -,727V associated with 
(j4^) ng N are defined inductively, for every p € {1, . . . , 27V}, by 

(4) f>= lim il(log||A* (A»_ 1 ...A%)\\), 

i=l 

Here, A P M denotes the pth exterior power of the matrix M, acting on the pth 
exterior power of M. 2N . One has 71 > . . . > 72jv- Moreover, the random matrices 
(A^)„ 6 n being symplectic, we have the symmetry property 72jv-i+i = — 7i, for 
every i£{l,...,N} (see 2 j p.89, Prop 3.2). 

To define Lyapunov exponents associated to the operator H(u>), we first introduce 
the sequence of transfer matrices associated to H(uj). Let E £ R and consider the 
second order differential system 

(5) H(uj)u = Eu. 

A function u = [u±, . . . , u^) : M. — > C N is called a solution of ((5|) whenever — u" + 
Vqu — Eu on R \ Z and each m satisfies the same boundary conditions as elements 
in D(H Ui ), that is 

(6) Vi € {1, . . . , N}, Vn E Z, u' l {n + ) = uj(n~) + Ul {n). 

Definition 3. If u is a solution of (0) then the transfer matrix Af n n _ ) _ 1 i (E) from 
n + to (n + 1) + is defined by the relation 

( u{{n + l)+) \ _ (») / u{n+) 

[i) y U '((n+l)+) J - A (n,n+l]K*)y 

Then for every n € Z, (J5) € Sp N (R) as*(u, u') is a solution of the first order 

Hamiltonian system associated to system (l5l). Thus the sequence ( Af 1 " 1 -.JE)) 

is a sequence of i.i.d. symplectic matrices and the Lyapunov exponents associated 
to it are by definition the Lyapunov exponents of H(uo). 

We have an explicit form for Af n ,,' i+1 ](-E)- To compute it we start by solving the 
free system ([5]) on (n,n + 1). Actually we only have to do it on (0, 1) due to the 
1-periodicity of Vo. Then the transfer matrix from n + to (n + I)~ is given by 

(8) A M (E) = exp ( Vq _° EIn j ; 
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We also set, for every N x N matrix Q the 2N x 2N matrix M(Q) given by 
M(Q) = ( 3? ? ) . Then, using the interface relation (jfjj), the transfer matrix 



Q In 

from (n + 1)~ to (n + 1) + is M (diag(cia;|™\ . . . , cjvuA^)). Thus we have 
(9) A^ n+1] {E) = M(diag(c 1 ^" ) , . . . ,c N J^)) A (0A) (E). 

The transfer matrix splits in a product of two factors. The first factor contains the 
random part of the transfer matrix and is independent of E. The second factor is 
deterministic and depends only on E. 

As the matrices A 1 ? ^ +1 ^(E) are i.i.d., we denote by /ie their common law. 

Definition 4. The Fiirstenberg group associated to H{oj) is the closed group gen- 
erated by the support of fiE ■' G{E) — < supp [Ie > • 

Because the transfer matrices are i.i.d. we have an internal description of G(E), 



(10) G(E) = < Afg^ (E) | € supp v >. 

2.2. Proof of theorem [TJ Using the criterion of Gol'dsheid and Margulis and 
Kotani's theory, theorem [1] reduces to the following proposition. 

Proposition 1. Let N = 3 and assume ci £ R and c\ and C3 being non-zero. 
There exists a discrete subset 5 o/l such that for every E £ R\<S , G{E) is Zariski 
dense in 



The rest of this section is devoted to the proof of this proposition. We begin by 
writing the explicit form of A^j, (E) . For this we have to compute the exponential 
in A( ^(i?) . We start by diagonalizing in orthonormal basis the symmetric matrix 
defining Vq : 







1 






(11) 


r 

V = 1 





1 


-( 




V 


1 









with 

(12) U = i I V2 -V2 I and U^ 1 =' C/. 




We now assume that i5 > \/2 and we will deal with the other cases later. By 

computing the successive powers of ( ^ n l one gets 

V Vq — EIn U I 



(13) A( ,i)(S) = I „ „ R a ,/3 



u \ „ /J/- 1 



[/ / a ' p ^ 7 V u 



T-i 



where a = y/E-1, (3 = V E - y/2, 7 = V E + ^2 and 



(14) R a ,nn 



1 


cos a 








— sin a 

a 










cos/3 








^ sin /3 










cos 7 










—a sin a 








COS OL 










—(3 sin/3 








cos/3 


V 








—7 sin 7 









\ 



- sin 7 



cos 7 / 
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To prove the Zariski denseness of G(E) in Sp 3 (R) we use the fact that Sp 3 (R) is a 
connected Lie group. If C\z(G(E)) denotes the Zariski closure of G(E) in Sp 3 (R) , 
we only have to show that 23(E), the Lie algebra of Glz{G(E)), is equal to sp 3 (R), 
the Lie algebra of Sp 3 (R) . We recall that 

(15) sp3(K) = I ( £ -* a ) ' a E -^ 3 ( K )' ^1 and h symmetric 

and sp 3 (R) is of dimension 21. Our strategy will be to exhibit a family of 21 
linearly independent elements in 23(E). Before starting this construction we prove 
the following lemma. 

Lemma 1. For a three-by-three matrix Q one has M(Q) £ C\z(G(E)) if and only 
if (I l)^s(E). 

Proof. Assume that ^ ° ° J £ g 3 (E). ThenM(Q) = exp(M(Q)-I) £ C\ Z (G(E)) 

because C\z(G(E)) is a Lie group. Conversely, if M(Q) £ C\z(G(E)) we consider 
the subgroup Gq := {M(nQ) = M(Q) n : n £ Z} of C\ Z (G(E)). It follows that 
M(xQ) £ CIz(Gq) for all x £ R. To see this, let p be a polynomial in 6 x 6 
variables such that p(A) = for all A £ Gq. Then the polynomial in one variable 
p(x) := p(M (xQ)) vanishes at every integer point, thus it vanishes identically on 
R. So, for every p(M(xQ)) = 0. Then by definition of the Zariski closure, 

for every x £ R, M(xQ) £ C1 Z (G Q ) C C\ Z (G(E)). By writing 

M(XQ) = I 2 N + X ( g q 

and by differentiating at the identity (at x = 0), by definition of a Lie algebra, 
Q o) e ^)- □ 

Proof of proposition^ for E > Step 1 : From ([9]), for and d)( ) in 

supp v we have 

(16) AfZWA^E)- 1 = Af(diag( Cl (^ 0) - ^ 0) ), . . .,c 3 (^ - J^))) 

is in G(E). Thus, 93 (E) being an algebra and therefore closed under linear combi- 
nations, using lemma[T]and the hypothesis that {x — y | x, y £ supp v} is spanning 

R 3 , we get that ( q ) ^ 03 (-^) f° r ever y diagonal matrix Q. If we assume that 

c 2 = then it is true only for Q of the form diag(ai, 0, 02), a±, a 2 G R. 

Step 2 : By step 1 and lemma[TJ M(Q) £ Clz(G(E)) for every diagonal matrix Q 

(or at least with a zero on the second diagonal term if c 2 = 0). In particular, 

(17) ^ (0> i)(^)=M(diag(ci4 0) ,---,c 3 4 0) ))" 1 Ao,i](^) G C1 Z (G(£)). 

Step 3 : By a general property in Lie groups, XMX^ 1 £ 23(E) for every M £ 
2s(E) and X £ G\z(G(E)). Thus, for every I £ Z, 

U \ i / \ , / t/- 1 
U J a 'P>~< \ U~ l QU y ^ J7 _1 

(18) A {0A) (E) 1 ( ° ° j^^-'ejJ^. 
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for any diagonal matrix Q. As U is orthogonal, 

U 
U 



(19) S3 (E) = sp 3 (R) & s 3 (E) 



So we are bring to prove that 83(E) = sp 3 (I 
I € Z and every diagonal matrix Q, by (fl8|) . 



(20) 



7? 





jj-^c/ 



33(E) 



u- 1 









= sp 3 (K) 



For this, we will use that for every 



K e fl 8 (B). 



Step 4 : One can choose Q = diag(v / 2, 0, \/2) (even for C2 = 0) to get 



U^QU 



For every 1 6 R, we set s(x) = sin(ar) and c(x) — cos(x). By step 3, for every I € Z 
we have 



( 



s(l"f)c(la) 
7 



-c(la)c(lf3) 
\ -c{la)c(li) 



s(la)c(lf3) 

OL 




-c(la)c(l(3) 








-e(Za)c(/ 7 ) 







s{la)s(lf3) 

a/3 
s(la) s (Z7) 



s(;a)c(;/j) 

s(Ze*)c(Z"y) 



s(Zo:)s(l/3) 



s(io:)s (Z7) 





0:7 














(iP)c(la) 







7 















\ 



is in g 3 (E). We can also choose Q = diag(2, 0, 2) (even if c 2 = 0) to get 

U- 1 QU = 




And so for every ! e Z, 

/ 2s(lu)c(lu) 
a 




2c 2 (la) 



Da (I) 





/3 

s(i7)c(i<3) 
7 







s(Wc(i 7 ) 



V 



c 2 (Z/3) 
c(J/3)c(J 7 ) 



2s" ! (ia) 




2c(la)s(la) 

OL 








s 2 (f/3) 

^~ 

s('/3)3(i7) 
/3 7 



s{10)c(10) 
/3 





jh, 



s('7)e(i/3) 
7 

s(i7)e(i7) 



c(//3)c(Z 7 ) 
c 2 (Z 7 ) 

is in 03(B). 

Step 5 : We prove that except for a discrete set of values of E, the matrices 
-Di(O), . . . , Di(7) are linearly independent. Indeed, if one computes the 8x8 de- 
terminant of the vectors 

/ — cos(la) cos(ip) \ 
cos(la) cos(/ 7 ) 



sin 


{la)sm(l{3) 


sin 


a/3 
[la) sin(/7) 


sin 


0:7 
la)cos(l/3) 


sin 


a 

la) cos(lf) 


sin 


l(3)°cos(la) 


sin 


l-y) cos(la) 



1 = 0,. 
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representing the D\(l) matrices, one gets : 

4096 sin 4 (a) sin 2 (/?) sin 2 (7) (cos a - cos /3) 4 (cos 2 (a) - cos 2 (/?)) x 
x (cos 2 (a) - cos 2 (7) ) (- sin 2 (2a) + cos 2 (f3) + cos 2 (7) 

(21) +2cos/3cos 7 (l - 2cos 2 (a))) 2 

which is a real analytic function of E on (s/2, +00) which does not identically 
vanish. Thus, this determinant vanishes only on a discrete set S\ C (v2) +00). 
Step 6 : Let E € (v2) +00) \ 6>i- By step 5 all matrices of the form 



(22) 



for (a,b,c,d,e, f,g,h) € K 8 are in 33(E). In particular we set Bo (respectively 
B\) the matrix of the form (|22p with a — 1 and the other parameters equal to 
(respectively b — 1 and the other parameters equal to 0). Then Bo £ 33(E), 
B x e 83(E) and 



f 


a 


6 





9 


h 


\ 


c 








9 










d 








h 













e 


/ 





— c 


-d 




e 








—a 










V / 








-6 








/ 



(23) 



and 



(24) 



Bo := [2?o,£i(0)] 



#a := [^,^(0)] = 



e 3 (^)- 



e 3 (^)- 



Step 7 : We prove that (B , B x , D 2 (0), . . . , D 2 (10)) is a family of 13 linearly 
independent matrices. For this we have to prove that a 13 x 13 determinant does not 
vanish identically. This determinant is given by the 11 vectors of the 13 different 
(and non-colinear) non-zeros elements in Z^O), Z^IQ) and the 2 vectors of 
the 13 corresponding coefficients in Bo and B\. Numerically one can verify that 
E = 1.6 S (V2, +00), E — 1.6 ^ Si and that for E = 1.6, the determinant has 
approximate value of —3507 =/= 0. Thus it does not identically vanish on (V2, +00). 
But this determinant being a real-analytic function of E on (\/2, +00) it is therefore 
vanishing only on a discrete set of values of E, S 2 C (V%, +00). 
Step 8 : Let S3 — Si U S2. Let E e (\/2, +00) \ S3. Looking at the zero 
coefficients in D\(l) and D 2 (l) one sees that the families (Di(0), . . . , D\(7)) and 
(Bo, Bi, D 2 (0), . . . ,D 2 (10)) lie in two orthogonal subspaces of sp 3 (R). Thus, they 
generate two orthogonal subspaces of dimension 8 and 13 and the direct sum of these 
spaces is still contained in 33(E). Thus, dim 33(E) > 21. But 33(E) C sp 3 (M) and 
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dim sp 3 (M) = 21, so Q3{E) = sp 3 (R). It endeed the proof for E > y/2 as by 
connectedness, Cl z (G(E)) = Sp 3 (R) for every E e (2, +oo) \ S 3 . □ 

Proof of proposition^ for E < y2- For E E (1, \/2) we have the same expres- 
sion for ^4(o,i) (E) as (fT3| but with changing a, (3, 7 into a — \JE — 1, [3 = y \f2 — E 

and 7 = V ' E + \/2. We also change in (fH| the cos(/3) into a cosh(/3), the term 
— /3sin(/3) into /3sinh(/3) and the term i sin(/3) into \ sinh(/3). Then we can follow 

the proof of the case E > y/2 which leads to a discrete set S4 c (1, V2) such that 
for every E £ (1, n/2) \ S 4 , C\ Z {G{E)) = Sp 3 ( R) . 

For E £ (-\/2, 1) wc set a = \/l - E, (3 = \/V2 - E and 7 = \J E + V2 and we 
do the same changes of cosinus and sinus into hyperbolic cosinus and hyperbolic 
sinus as in the case E £ (1,V2), for those involving a and (3. Then following the 
proof of the case we get the existence of a discrete set S5 C (— v2> 1) such 

that for every £ e (—v/2, 1) \<S 5 , C1 Z (G(£)) = Sp 3 (R) . 

For E G (-00, — v/2) we set a = VT^E, (3 = y/V^-E and 7 = V-S - \/2 and 
we change all the sinus and cosinus into hyperbolic sinus and cosinus in (jT4j) . It 
leads to a discrete set S$ c (—00, — v2) such that for every _B £ (—00, — V2) \ 56 , 
Clz(G(E)) = Sp a (R). 

Finally if we set S7 = {— -s/2, 1, V^}, by setting S = S3 U . . . U 1S7 we have a discrete 
set such that for every £ £ l\5, G(E) is Zariski dense in Sp 3 (M) , which proves 
proposition [TJ □ 

To finish the proof of theorem [1] it just remains to check that Kotani's theory can 
be applied for H (uj) which will be done in Section 13.2.21 Then from the non- 
vanishing of the Lyapunov exponents outside of iS, which is of Lebesgue measure 
0, we deduce the absence of absolutely continuous spectrum of H(iv). And thus 
theorem Q] is proved. 

In our study of the Lyapunov exponents it would be also interesting to look at 
what happens at the values of E in S. We do not know exactly if the Lyapunov 
exponents are all vanishing at these energies or if only some of them vanish or even 
if they are all positive but not distinct. To handle directly with the computation 
of the Lyapunov exponents at the values of E in S is much more difficult than in 
the case of scalar- valued operators. For the scalar- valued Anderson model in [B], 
the Lyapunov exponent vanishes at the critical energies, which is clear because at 
these values, G(E) is compact. Here it is no longer the case and except doing some 
numerical attempts to understand the situation, no rigorous proof has been found 
yet. We think to be in the presence of energies at which some Lyapunov exponents 
vanish while others are not. The case were all the Lyapunov exponents vanish does 
not seem to happen. 



3. Regularity results on Lyapunov exponents and the IDS 

3.1. Holder continuity of the Lyapunov exponents. In this section we prove 
theorem [3l We will see how to use general results of [4] for the operator H(uj). 
We can deduce theorem [3] from theorem 1, p. 885 in [4j once we have proved the 
following estimates on the transfer matrices. 
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Proposition 2. Let I be a compact interval in R. There exist C\ > 0, C2 > 
independent of n, oj( n \ E such that for every p G {1, . . . , N}, 

(25) || A p A^ +1] (E)\\ 2 < cMpCi + p\E\ + p) < C 2 . 

And there exists C3 > independent of n, u'™', E such that for every E, E' S I and 
every p € {1, . . . , N}, 

(26) ||A^^ +1] (£)-A^£ +1] (^)|| <C 3 \E-E'\. 

Proof. First we recall that -A(o,i) is obtained by solving the free system H{lo)u = 
Eu on (0, 1). This system is of the same form as the one from which we deduce the 
transfer matrices in [4]. Thus, the estimates (|25p and (|26|) which were proved for the 
transfer matrices in [3] are still valid for ^4(o,i) (E) and thus j4(o,i) (E) verifiy (f2"5)l and 

(f26|) . Then suppi/ being bounded, | |M(diag(ciuj™\ . . . , cnlo^))\ \ will be bounded 
uniformly in n, w' n ) and also E because it does not depend on E. Let C > 
be independent of n, a/™' and £7 such that \\M(diag(ciUj[ n \ . . . ,CNto^))\\ < C. 
Then || A p M(diag(ci^ n) ,.. .,c N ^P))\\ < C p . Using the fact that || A p (MiV)|| < 
1 1 A p M 1 1 1 1 A p N\ I for every matrices and the relation |9]), we finally obtain (|25|) and 
®for^ +1] (£J). □ 

Applying theorem 1 in [4] and proposition [5] we have proved theorem [3] Then 
applying theorem [3] on every compact interval I C R \ S where S is obtained in 
theorem [TJ we get theorem HI 

3.2. Holder continuity of the IDS. 

3.2.1. Existence of the IDS. Once again we will follow the method used in [4] to 
prove this time the existence of the integrated density of states associated to H(ui) 
and its Holder continuity. The proof of the existence of the limit ([3|) is based 
upon the fact that the one-parameter semi-group (e~ tR{D) ^)t>o admits an integral 
kernel in L 2 (D 2 ) ® _Mn(C). This kernel is coming from a Feynman-Kac formula. 
There already exists such formula for scalar-valued point interactions operators 
as presented in [?]■ Adapting the Borel measure representation method of [7], 
using Lie- Trotter formula as it is done in [4] and noticing that the time-ordered 
exponential in 4|] becomes now a usual exponential, one gets 

(27) V/ e L 2 (R) ® C", Vx € R, e' tH ^f(x) = f K t (x, y)f(y)dy 

JR 

with 

(28) Vt > 0, Va, yel, K t (x, y) = / e~ tv °d^ x , y ^(w) 



where for every to G fixed, ^. x .y.Lo is a finite measure on the space of the 

continuous paths w on [0, t] such that w(0) = x and w(l) = y. Then to deduce the 
kernel of e~ tHi we introduce Td(w), the time of first exit from D of the path 

W € \N x ,y, 

(29) Td(w) = inf{t > I w(f) ^ £>}. 
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Then we have (see [12] ) 

(30) VfeL 2 (D)®C N , VarGD, e- tHiD) ^f{x) = / K[ D \x,y)f{y)Ay 
with 



(31) Vt > 0, Var.j/ G D, ^(^y) = 6 " / ^^(w^^d^^w). 

V 27TT J 

As we can see, for every t > 0, if t is in L 2 (D 2 )® .Mn(C). Thus all the discussion 
made in Section 2.3 of [I] applies here to get the existence of the IDS for every 
E G M and its realization as the distribution function of a measure n called the 
density of states. 

3.2.2. Kotani 's theory and proof of theorem^ We start by adapting Kotani's theory 
of [2] to our setting. According to the presentation made in [3] we actually only 
have to prove that theorem 2.1 (a) of [2] is true for the operator H(ui). We fix 
lu en. Let C + ={z£ C, Imz > 0} and C_={z£ C, Imz < 0}. For E G C+UC_ 
we set 



\f(x)\ 2 dx < +oo 



and 



f N r- 

J+(E) = { f 6 0J)(ffJ | = £/ and / 

J_(£) = |/e®D(fl Ul ) | H(u)f = Ef and J \f(x)\ 2 dx < +ooj . 
Proposition 3. We /iaue : dim J+(E) = dim J-(E) = N. 

Proof. Each fl a)j is in the limit point case of singular Weyl theory (see pQ) which 
ensure the existence of a unique solution f\ t i S D(H m ) of = -E 1 / such that 

Jq°° \ fi,i( x )\ 2 dx < +oo and a unique solution /2.i G D(H Wi ) of = J?/ such 

that f° oD \f 2 , i (x)\ 2 dx < +oo. Then if we introduce H+ = ) n i, 2 (o,+oo), 
dimKer(_ff+ — E) = 1. With H~ = -Ho)Jz3(H„.)ni 2 (-oo,o)) we also have dim Ker (JI^ - 
i?) = 1. Then 



dim Ker ^® ff+ - E^j = dim Ker ^® H ~ - Ej 



N. 



As Vb is bounded, from deficiency index theory we get 

dim J+{E) = dim Ker ^® H+ + V - Ej = dim Ker ^® H+ - E^j = N 

and the same goes for J -(E). □ 

From this proposition we get as in [T3], Corollary 2.2, the following result. 

Proposition 4. Let E G C+ U C_ and lu G fl. Then there exists a unique function 
x i — ► F + (x,E) with values in A^n(C) (respectively x i— > F-(x,E)) satisfying 

/>oc 

H(ui)F + = EF + , F+(Q,E) = I, and / £T)|| 2 dx < +oo, 
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respectively 



H(u>)F- = EF-, F-(Q,E) = I, and 



\\F-{x, E)\\ 2 dx < +00 



—00 



This proposition is the starting point of all the theory on the Floquet exponent w 
and the M±-functions presented in pH] , Thus Kotani's theory on the absolutely 
continuous spectrum apply for H{uS) and it finishes the proof of theorem [T] We also 
recall that [14] only considers R-ergodic operators whereas H(ui) is Z-ergodic. To 
avoid this difficulty we can refer to the suspension procedure developed by Kirsch 
in This procedure allows us to construct from H(u>) an operator H(uj) define 
on a bigger probability space which is M-ergodic. H(lo) is also constructed in a way 
such that its IDS and Lyapunov exponents exist if and only if those of H(u>) exist 
and in this case they are equal for both operators. 

Also, from the properties of the Floquet exponent w combined with previous results 
of Kotani (see [13]) one can repeat the discussions of Sections 4.1 and 4.2 in [1] to 
prove the following Thouless formula for H(u). 

Proposition 5. For Lebesgue- almost every E el we have 



where a is a real number independent of E and n is the density of states of H(lu). 
Moreover, if I Cl is an interval on which E 1— > (71 + . . . + r jM)(E) is continuous 
then ([Hp holds for every Eel. 

Using this Thouless formula, theorem [3] and properties of the Hilbert transform 
(see |15j). we can obtain the proof of theorem [5] exactly in the same way as it was 
done at theorem 4 of [4], Section 4.3. Then applying theorem [5] on any compact 
interval / cK\S where S is given in theorem [TJ we prove theorem [Fj] 

As we can see, the methods to prove regularity of Lyapunov exponents and regu- 
larity of the IDS for the point interactions model H(u) are completely similar to 
those for matrix- valued Anderson models. The main differences are to be found 
in the proof of the Zariski denseness of the Fiirstenberg group. Indeed, in [3] we 
proved Zariski denseness of the Fiirstenberg group of an Anderson operator acting 
on L 2 (M) ® C 2 using algebraic technics different than those used in the proof of 
theorem [1] in the present paper. 
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